Quasispecies evolution on a fitness landscape with a fluctuating peak 
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A quasispecies evolving on a fitness landscape with a single peak of fluctuating height is studied. In 
the approximation that back mutations can be ignored, the rate equations can be solved analytically. 
It is shown that the error threshold on this class of dynamic landscapes is defined by the time average 
of the selection pressure. In the case of a periodically fluctuating fitness peak we also study the 
phase-shift and response amplitude of the previously documented low-pass filter effect. The special 
case of a small harmonic fluctuation is treated analytically. 
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I. INTRODUCTION 



II. QUASISPECIES IN DYNAMIC 
ENVIRONMENTS 
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Ever since Eigen's work on replicating molecules in 
1971 [Q, the quasispecies concept has proven to be a 
very fruitful way of modeling the fundamental behavior 
of evolution. A quasispecies is an equilibrium distribu- 
tion of closely related gene sequences, localized around 
one or a few sequences of high fitness. The kinetics of 
these simple systems has been studied in great detail, 
and the formulation has allowed many of the techniques 
of statistical physics to be applied to replicator and evo- 
lutionary systems, see for instance [Q-||. 

The appearance in these models of an error-threshold 
(or error-catastrophe) as an upper bound on the muta- 
tion rate above which no effective selection can occur, has 
important implications for biological systems. In partic- 
ular it places limits on the maintainable amounts of ge- 
netic information jl]|2| which puts restrictions on theories 
for the origins of life. 

Until now studies of quasispecies have mainly focused 
on static fitness landscapes. However, many organisms in 
nature live in a quickly changing environment. In this pa- 
per we will study how a population responds to changes 
in the fitness landscape. More precisely we will study 
the population dynamics on a fluctuating single peaked 
fitness landscape. Since the full theory turns out to be 
impossible to solve analytically, we introduce a simple ap- 
proximation that makes the rate equations analytically 
tractable. The expression for the error threshold is then 
obtained from the expression in the static case by replac- 
ing the height of the fitness peak by the time average of 
the height of the fluctuating peak. We also study how 
the phase-shift between fitness oscillations and popula- 
tion dynamics depends on the frequency in the case of a 
small harmonic fluctuation. 



A quasispecies consists of a population of self- 
replicating genomes, where each genome is represented 
by a sequence of bases sjt, [s\S2 ■ ■ • s v ). We assume that 
the bases are binary, s k € {1,0} and that all sequences 
have equal length v. Every genome is then given by a 
binary string (011001 • ■ ■), which also can be represented 
by an integer k = £\ Sj2 j (0 < k < T). 

To describe how mutations affect a population we de- 
fine Wf. as the probability that replication of genome I 
gives genome k as offspring. We only consider point mu- 
tations, which conserve the genome length v. 

We assume that the point mutation rate \i = 1 — q 
(where q is the copying accuracy per base) is constant 
in time and independent of the position in the genome. 
We can then write down an explicit expression for W\ in 
terms of the copying fidelity: 



(1) 



where hki is the Hamming distance between genomes k 
and I. The Hamming distance hki is defined as the num- 
ber of positions where the genomes k and I differ. 

The equations that describe the dynamics of the pop- 
ulation take a relatively simple form. Let Xk denote the 
relative concentration and Ak(t) the time-dependent fit- 
ness of genome k. We then obtain the following rate 
equations: 



x k {t) = ^W{Ai{t)xi{t) - e{t)x k {t) 



(2) 



where e(t) = Ai(t)xi{t), and the dot denotes a time 
derivative. The second term ensures the total normal- 
ization of the population (%2i x l(t) = 1) so t na t x k(t) 
describe relative concentrations. 
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In the classical theory introduced by Eigen and cowork- 
ers [00.1^], the fitness landscape is static. The rate equa- 
tions (p ) can then be solved analytically by introducing a 
change of coordinates that makes them linear, and then 
solving the eigenvalue system for the matrix W^Ai. The 
equilibrium distribution is given by the eigenvector cor- 
responding to the largest eigenvalue. 

If the fitness landscape is time-dependent, this method 
cannot be applied. A time-ordering problem occurs when 
we define exponentials of time-dependent matrices, since 
in general the matrix W^Aiif) does not commute with 
itself at different points in time. Later in this paper we 
make a simple approximation that makes the rate equa- 
tions one-dimensional; time-ordering is then no longer 
necessary. 

Much of the work on quasispecies has focused on fitness 
landscapes with one gene sequence (the master sequence) 
with superior fitness, a, compared to all other sequences. 
These are viewed as a background with fitness 1. These 
landscapes are referred to as single peaked landscapes. 
The master sequence is denoted xq. In this paper we fo- 
cus on single peaked landscapes where the height of the 
fitness peak is time-dependent. The fitness landscape is 
then given by 



A k (t) 



a{t) if k = Q 
1 otherwise 



(3) 



This class of time-dependent landscapes was studied by 
Wilke and co-workers • They investigated the be- 

havior of a periodically fluctuating single peak landscape 
by numerically integrating the dynamics to find the limit 
cycle of the concentrations for a full period. 

Fig. |l| shows how the concentration of the master se- 
quence responds to a sudden, sharp jump in its fitness. 
When the fitness changes it takes some time for the pop- 
ulation to reach the new equilibrium. It is this delay 
that causes a phase shift between a periodically changing 
fitness function and the response in the concentrations. 
The relaxation time of the population to the appropri- 
ate equilibrium distribution depends on both the fitness 
values of the landscape and the mutation rate. For ex- 
tremely slow and smooth changes in the fitness the popu- 
lation will effectively reach equilibrium at every point in 
time. Thus the continued existence of a quasispecies will 
depend on the local dynamics of the landscape. When 
the landscape changes quickly the population will fail 
to follow the changes adequately and thus responds to 
the landscape dynamics in a way that is typical of a low 
pass filter. The following section examines the fluctuat- 
ing single peak landscape in some detail. In particular, 
we introduce an approximation that lets us find an ana- 
lytic form for the relaxation time of the population, and 
the phase lag it introduces in a periodic landscape. 
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10 11 12 13 14 15 

FIG. 1. The concentration of the master sequence when 
the fitness peak makes a sudden jump. The fitness is given 
by a — 10 when t £ [10, 12], a = 5 otherwise. The genomic 
copying fidelity is given by Q — 0.7. 



III. APPROXIMATE QUASISPECIES DYNAMICS 

We now introduce a simple approximation of the model 
presented above. In this approximation we can solve the 
rate equations and find an expression for the concentra- 
tion of the master sequence xo(t). In the limit of long 
chain-length (v ^> 1) we can neglect back- mutations from 
the background to the master sequence. This gives a sim- 
plified one-dimensional version of the rate equation of the 
following form: 



&o(t) = Qa(t)x a (t) - e{t)x (t) 



(4) 



where Q = q v is the copying fidelity of the whole genome 
and e(t) = (a(t) - l)x (t) + 1. 

Fig. p] compares the concentration of the master se- 
quence calculated by solving approximation \\ and by nu- 
merically integrating the full rate equation [A The figure 
shows that the approximation is quite accurate. 
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FIG. 2. The concentration of the master sequence calcu- 
lated using the approximation ^ (dashed) and by numerically 
solving the rate equations ^ (solid). The fitness is given by 
a = 10 + 5sin(4t). The genomic copying fidelity is given by 
Q — 0.7 and the genome length v = 25. 

Since this equation is one-dimensional there is no time- 
ordering problem and it can be solved analytically for 
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non-periodic peak fluctuations. Equation Q can be trans- 
formed to a linear form by introducing a new variable 



y(t) = 1 - (Qa(t) - l)y(t) 



(5) 



which can be solved. Substituting back gives the concen- 
tration of the master sequence 



x (t) = x - 



I rs 



(6) 



Since we are only interested in the long time behavior 
of the system we can ignore transients carrying memory 

from initial values. Assuming ^ <7 ^ 1 ^ ds ;g> l gives 



x Q (t) = 



Q 



l + (l-Q)foe-J> {Qa{u) - 1)du d S 



(7) 



This is a generalization of the static expression for the 
asymptotic concentration: 



Qa- 1 

a - 1 



(8) 



On a static single peaked fitness landscape there is a 
phase transition in the concentration distribution when 
the copying fidelity decreases below a critical value [§,|l3| . 
At high mutation rate the selective advantage of the 
master sequence due to its superior fitness is no longer 
strong enough for the gene sequences to be localized 
in sequence space. Instead they diffuse over the entire 
sequence space, and the distribution becomes approxi- 
mately uniform. This is generally referred to as the error 
catastrophe or error threshold and is one of the main im- 
plications of the original quasispecies model. By making 
the same approximation as above, i.e. assuming no back- 
mutations onto the master sequence, the static landscape 
error threshold can be shown to occur when Q = -. In 
other words, the transition occurs when the selective ad- 
vantage of the master-sequence no longer is able to com- 
pensate for the loss of offspring due to mutations. This 
can also be seen from Eq. |S| which defines the stationary 
distribution of the master sequence in the static case. 

One has to be careful when discussing the error thresh- 
old on a fluctuating peak. The fitness can, for example, 
slowly move from being strong enough to localize the 
population around the peak, to beibg so weak that the 
population delocalizcs, and then back again. If we how- 
ever consider an average over a time scale much longer 
than the fluctuation time of the fitness peak, a sensible 
definition of the error threshold can be made based on the 
average concentration of the master sequence. The time 
average of the concentrations can be found by rewriting 
equation [| as differentials 



dxo 
x 



(Qa(t) - 1 - (a(t) - l)x (t))dt 



(9) 



The concentration of the master-sequence is positive. 
The left hand side of Eq. ^ is therefore positive and 
the last term in the integral, — (a(t) — l)xo(t), on the 
left hand side is negative. This implies that for xo(t) 
to be positive as time goes to infinity, we must assume 
/ (Qa(t) — 1) dt > 0. The fluctuating time dependent 
equivalent to the static error threshold is therefore given 
by 



Qfc 



1 

Wt 



(10) 



This shows that the error threshold on a fluctuating fit- 
ness peak is determined by the time average of the fitness, 
if the fluctuations are fast compared to the response time 
of the population. 

Eq. fj] indicates that the response time of the system 
is approximately given by (Qa(t) — 1) , i.e. the rela- 
tive growth of the mastersequence compared to the back- 
ground. For the time average mentioned above to be an 
interesting parameter the fluctuations of the fitness peak 
must therefore be faster than this response time; only 
for this kind of environmental dynamics is it sensible to 
talk in terms of the average concentration of the master- 
sequence. Thus if the fluctuations occur on a time-scale 
faster than the response-time of the quasispecies, then 
the error-threshold is defined by Eq. |l^. For extremely 
slow changes the system will effectively be in equilibrium 
around the current value of the fitness. For slightly faster 
changes the response of the population will lag somewhat 
behind the changes in selective environment. In these 
cases it is more interesting to study the minimal concen- 
tration of the master sequence, which occurs when the 
fitness peak has a minimum (as we shall see later the 
phase-shift decreases when the fluctuation frequency de- 
creases). 

When the full replicator equations for a rapidly fluctu- 
ating peak are numerically integrated, the time-averaged 
quasispecies distribution displays an error catastrophe at 
high error rates fi = 1 — q. In figure || the fitness peak 
fluctuates periodically with <j(t) = 10 + 5sin(i). The av- 
erage fitness is given by (a)t = 10 and the genome length 
v = 25 and thus Eq. |l0] predicts the error-threshold to 
occur at fi = 0.088, which agrees with the value found 
by numerically integrating the equations of motion di- 
rectly. The analysis in this section demonstrates that 
by making the error tail approximation and reducing the 
dynamics to one-dimensional form, an analytic form ex- 
ists for the error-threshold on fast moving landscapes. 
This one-dimensional formulation removes the need to 
time-order the changes in selective advantage of the land- 
scape. This allows the integrals for the time history of 
the master-sequence concentration to be solved explicitly 
in equation 0. 
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FIG. 3. The time-averaged quasispecies distribution is 
shown as a function of the error rate /x = q — 1. The figure 
shows the numerical solution to the full rate equations. The 
fitness peak is defined as a(t) = 10 + 5sin(t) and the genome 
length is v = 25. The error threshold is located at /j, m 0.085, 
corresponding to Qf c — 0.109 which can be compared to the 
approximate value Qf c — 0.1 predicted by Eq. nO. 



IV. PHASE-SHIFTS ON PERIODIC 
LANDSCAPES 



To study how the master sequence responds to changes 
in the height of the fitness peak it is convenient to assume 
that the fluctuations are periodic. It then makes sense 
to speak of the amplitude of the oscillations in concen- 
tration and of the phase-shift between the concentration 
and the fitness. It is intuitively clear that when the fit- 
ness peak is oscillating slowly (compared to the response 
time (Q<r(t) — l)^ 1 ) there will be a very small phase-shift; 
the population will have time to reach an equilibrium 
about every value of <r(t). The amplitude of changes in 
the master-sequence concentration will, for the same rea- 
son, be as large as possible. This result, together with 
the time- aver aging effect found in the previous section, 
indicates that the population responds to the driving of 
the environment with a low pass filter effect. In one- 
dimensional population genetic models this phenomenon 
has been noted for some time [O Wilke et al. |Tl| l 
demonstrated via simulations that the same filtering oc- 
curred to quasispecies evolution on a periodically fluctu- 
ating single peak. Noting that the maxima and minima 
in concentration occurs when xq = 0, we can find a rela- 
tion between the phase-shift (between the concentration 
and fitness fluctuations), and the amplitude of the fitness 
fluctuations. Let t xmax be the time when the concentra- 
tion has a maximum. Similarly the fitness is at a max- 



imum at time t, 
two is 5 = t x 



a max • 



t„ 



Thus the phase-shift between the 
lax . From equation ^ the condition 
for the maximum value of xq during a full cycle can be 
derived 



a max + S)-l 
max(x (*)) = —77 --r — - 



time (t) 

FIG. 4. Thhe response in concentration of the master se- 
quence (solid line) as the fitness peak oscillates according to 
o(t) = 10 + sin(4t). The genomic copying fidelity is Q = 0.7. 
The dashed line shows cr(t), scaled to fit in the plot. Note the 
phase-shift between the fitness function and the concentration 
response. 

In general there is no closed analytic expression for this 
phase-shift (6) , or the response amplitude of the master- 
sequence concentration. When the fluctuations of the 
fitness peak is a small harmonic oscillation equation |ll| 
becomes analytically tractable. For such fluctuations 

y{t) = 1 - (Qa(t) - l)y(t) (12) 
a(t) = a + e sin(wi) (13) 

From equation (7] it is reasonable to assume the solution 
to be of the form y(t) = q^_ x + u(t), where u(t) is small 
compared to the average. Ignoring higher order terms 
equation || can be written in terms of the perturbation 
u(t) as 



u(t) = (1 - Qa)u(t) 



! sin(wt) 



Qa-1 

This differential equation can be solved to obtain 
eQ 



(Qa-l)V(Q<r-l) 2 

LU 



(14) 

(15) 
(16) 



u(t) = —77^3 — — 7== sm(ujt — 6) 

tan(<5) 

Qcr - 1 

In eq. [l5| and |l6| transients have been ignored since they 
decay exponentially as e~( < 2 <T-1 - ) *. Thus the frequency of 
the oscillations is normalized by the (average) response 
rate of the population Qa — 1. 



(11) 
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FIG. 5. The phase shift as a function of the period T = — . 
The dashed line is a prediction using Eq. and the solid is 
derived by numerically solving the reate equations Param- 
eters used are a = 10 + sin(urf), Q = 0.7 and v = 25. 

Substituting this back into the expression for xo(t) 
gives 

x 

(ff-i)V(Qff-i) 2 +" 2 

where x — _ . 

(7 — 1 

The characteristic behavior of a low pass filter is clearly 
shown in equation |l^ and |l7]. As the frequency of the 
fluctuations increases, the amplitude of the concentra- 
tion response decreases and the phase shift converges to 
Figure |^ shows how a population responds to har- 
monic oscillations of the fitness peak. The phase-shift 
makes the concentration of the mastersequence reach its 
maximum when the actual fitness has already decreased 
below maximum. 



V. CONCLUSIONS 

In this paper we have shown that the time dynam- 
ics of a quasispecies on a fluctuating peak can be stud- 
ied under the standard no back-mutation approximation. 
The general time ordering problem stemming from a time 
dependent landscape disappears since the rate equation 
becomes one-dimensional. We show that the time depen- 
dent equivalent to the static error threshold is determined 
by the time average of the fluctuations of the fitness peak. 
An expression for the typical response time for a popu- 
lation is given in terms of copying fidelity and selection 
pressure. We also show that for small periodic fluctua- 
tions the time dynamics of the population has a phase 
shift and a low pass filter amplitude response. Analytic 
expressions for the phase shift and the amplitude are de- 
rived in the special case of small harmonically oscillating 
fluctuations. 
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